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Abstract 

Previously,  expressions  governing  the  temporal  evolution  of  linear  perturbations  to  an  isolated, 
planar,  two-dimensional  shock  front  in  an  inviscid  fluid  medium  with  an  arbitrary  equation  of 
state  were  derived  using  a  methodology  based  on  Riemann  invariants  and  Laplace  transforms 
[J.W.  Bates,  Phys.  Rev.  E  69,  056313  (2004)].  An  overlooked  yet  immediate  consequence  of  this 
theory  is  that  the  stability  limits  of  shocks  can  be  readily  determined  from  a  inspection  of  the 
poles  of  the  transformed  ripple  amplitude.  Here,  it  is  shown  that  two  classes  of  instabilities  exist 
for  isolated  planar  shock  waves:  one  in  which  perturbations  grow  exponentially  in  time,  and  the 
other  in  which  disturbances  are  stationary.  These  results  agree  with  those  derived  by  D’yakov 
and  Kontorovich  (by  more  arduous  and  somewhat  ambiguous  means),  and  serve  as  an  important 
addendum  to  our  earlier  analysis. 
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I.  INTRODUCTION 


The  stability  of  shock  waves  is  an  issue  of  fundamental  importance  in  compressible  fluid 
dynamics,  with  far-reaching  implications  in  fields  ranging  from  astrophysics  [1]  to  inertial 
confinement  fusion  [2],  Under  most  circumstances,  a  small  sinusoidal  disturbance  to  an 
otherwise  planar  shock  front  in  a  homogeneous  fluid  will  oscillate  with  an  amplitude  that 
decays  asymptotically  according  to  an  inverse  power  law  ( e.g .,  as  t-3/2,  where  t  denotes  time 
[3-5]).  Even  in  the  absence  of  viscous  damping  effects,  all  linear  disturbances  to  the  front 
will  eventually  die  away,  and  the  shock  will  regain  its  planarity.  Roberts  [6]  was  the  first  to 
show  that  perturbed  shocks  in  fluids  obeying  an  ideal-gas  equation  of  state  (EOS)  always 
display  this  asymptotic  behavior,  and  are  thus  unconditionally  stable.  For  a  fluid  with  a 
non-ideal  EOS,  however,  shock  stability  is  not  guaranteed,  and  under  certain  conditions, 
small  perturbations  actually  may  amplify  rather  than  evanesce,  ultimately  leading  to  a 
reorganization  of  the  steady-state  flow  [7-10]. 

The  initial  theoretical  investigation  of  this  subject  was  performed  over  50  years  ago  by 
D’yakov  [11],  who  used  a  normal- mode  analysis  to  determine  conditions  in  which  “corruga¬ 
tions”  on  an  isolated  shock  traveling  through  an  arbitrary  quiescent  fluid  at  constant  speed 
D  would  grow  exponentially  in  the  linear  approximation.  The  use  of  the  term  “isolated” 
here  refers  to  an  idealized  situation  in  which  a  shock  is  decoupled  completely  from  its  driving 
mechanism  ( e.g .,  a  moving  piston,  or  an  ablation  surface).  Although  derived  in  a  somewhat 
cryptic  fashion  with  reliance  of  numerous  variable  transformations,  D’yakov’s  results  proved 
to  be  seminal  because  they  stated  for  the  first  time  shock  instability  criteria  in  terms  of  sim¬ 
ple  physical  quantities  that  were  useful  to  experimentalists.  Using  the  so-called  “D’yakov 
parameter” 

V  \dpj  H 

whose  derivative  is  taken  along  the  Hugoniot  curve  in  the  pressure-density  (p,  p )  plane,  and 
evaluated  at  the  downstream  state,  these  criteria  can  be  summarized  as 

h  <  —  1  or  h  >  1  +  2 Mi .  (2) 

Here,  Mi  =  (D  —  U)/c\  is  a  downstream  Mach  number  satisfying  0  <  Mi  <  1,  where  U 
is  the  shocked  fluid  velocity  in  the  laboratory  frame,  and  Ci  is  the  associated  sound  speed. 
The  quantity  r)  =  Pi/po  >  1  in  Eq.  (1)  is  the  compression  ratio  across  the  shock,  where  the 
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subscripts  “0”  and  “1”  denote  upstream  and  downstream  states,  respectively. 

In  the  1970s,  Swan  and  Fowlcs  [12,  13]  revisited  the  shock  stability  problem  in  an  effort 
to  clarify  several  issues  that  arose  in  D’yakov’s  original  analysis.  These  authors  carefully 
repeated  the  laborious  normal-mode  calculation,  and  found  that  despite  the  ambiguities 
in  his  derivation,  D’yakov’s  conclusions  were  nevertheless  correct.  Furthermore,  Swan  and 
Fowles  showed  that  the  hirelings  of  a  separate  investigation  due  to  Erpenbeck  [14]  —  whose 
mathematical  approach  was  different  and  whose  results  were  cast  in  terms  of  thermodynamic 
quantities  —  were  equivalent  to  the  inequalities  appearing  in  Eq.  (2).  Using  Erpenbeck’s 
formulation  of  the  instability  criteria,  Gardner  [8]  was  able  to  show  that  satisfaction  of 
the  condition  1  +  2M\  <  h  leads  to  the  splitting  of  a  shock  into  two  counter-propagating 
waves  —  although  apparently  this  phenomena  has  never  been  observed  experimentally.  The 
other  instability  criterion,  h  <  —1,  is  known  to  correspond  to  the  break-up  of  a  front  into 
two  waves  traveling  in  the  same  direction  [7],  and  has  been  documented  for  some  time  in 
experiments  with  shocked  materials  undergoing  phase  transformations,  or  yielding  at  the 
elastic  limit  [15]. 

In  addition  to  temporally  exponentiating  solutions,  D’yakov’s  analysis  (which  contained 
an  error  that  was  later  corrected  by  Kontorovich  [16])  also  identified  a  regime  in  which  the 
amplitude  of  shock  perturbations  remains  stationary.  This  phenomenon  is  known  as  the 
“D’yakov-Kontorovich”  (DK)  [17,  18]  or  “acoustic  emission”  [19,  20]  instability,  although 
such  terminology  is  arguably  misleading  since  the  growth  rate  is  zero  for  all  perturbative 
modes,  and  no  true  instability  takes  place  —  at  least  not  in  the  conventional  sense.  Once 
created,  small  disturbances  on  the  front  persist  indefinitely  —  without  amplifying  or  damping 
—  and  continuously  radiate  sound  and  entropy  vortex  waves  into  the  downstream  flow  behind 
the  shock.  The  DK  instability  can  occur  when  h  satisfies 


hc 


1  -  M*(  1  +  r,) 
1  -  Mf(l  -  r) ) 


(3) 


Fowles  [21]  has  shown  that  this  same  expression  can  be  derived  by  determining  the  conditions 
for  the  reflection  coefficient  of  acoustic  waves  incident  on  a  shock  front  from  the  downstream 
direction  to  become  infinite.  Note  that  the  parameter  hc  in  Eq.  (3)  always  lies  between 
the  limits  —1  and  1,  therefore  extending  the  range  of  values  of  h  for  which  “unstable” 
shock  behavior  may  occur.  For  an  ideal  gas,  hc  =  (1  —  2Mo)-1  and  h  =  —1/Mg  (where 
M0  =  D/cq  >  1  is  the  Mach  number  of  the  shock  in  the  upstream  fluid,  and  cq  is  the  speed 
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of  sound  there),  and  so  the  inequality  in  Eq.  (3)  can  never  be  satisfied.  There  is  evidence 
to  suggest,  though,  that  planar  shocks  in  other  substances  may  experience  this  type  of 
instability.  Examples  include  van  der  Waals  fluids  [18,  20],  ionizing  [22,  23]  and  dissociating 
[24,  25]  gases,  as  well  as  certain  metals  [26]. 

In  this  paper,  we  present  an  alternative  derivation  of  the  shock  instability  criteria  appear¬ 
ing  in  Eqs.  (2)  and  (3)  that  is  more  tractable  than  the  traditional  normal-mode  approach. 
Although  our  conclusion  is  not  new,  the  means  that  we  employ  to  reach  it  is.  Our  analysis 
relies  on  a  methodology  based  on  Riemann  invariants  and  Laplace  transforms  that  was  pi¬ 
oneered  by  Roberts  [6],  and  extended  recently  [27]  to  account  for  a  fluid  with  an  arbitrary 
EOS.  The  salient  aspects  of  this  underlying  theory  are  summarized  in  Sec.  II,  although  we 
should  remark  that  in  the  interest  of  brevity,  most  of  the  mathematical  details  have  been 
omitted  and  only  expressions  directly  relevant  to  our  purposes  here  have  been  reproduced; 
consequently,  the  reader  is  encouraged  to  consult  Ref.  [27]  in  concert  with  what  follows. 
One  of  the  principal  results  of  our  earlier  investigation  was  the  derivation  of  an  equation 
for  the  Laplace-transform  of  the  shock-ripple  amplitude,  and  we  show  in  Sec.  Ill  how  the 
shock  stability  limits  can  be  gleaned  from  an  inspection  of  the  poles  of  this  expression  us¬ 
ing  a  simple  graphical  technique.  The  validity  of  these  limits  is  verified  also  in  Sec.  Ill  by 
comparison  with  numerical  solutions  of  the  Volterra  equation  that  the  ripple  amplitude  is 
known  to  satisfy.  Finally,  in  Sec.  IV  we  summarize  the  conclusions  of  this  study. 


II.  A  BRIEF  SUMMARY  OF  OUR  EARLIER  WORK 

Consider  an  isolated,  planar,  step  shock  as  described  in  Sec.  I  that  propagates  in  the 
negative  x- direct  ion,  and  has  a  front  of  infinite  extent  that  lies  parallel  to  the  y-axis.  In  the 
frame  of  reference  in  which  the  downstream  fluid  is  at  rest,  the  position  of  the  unperturbed 
shock  as  a  function  of  time  is  given  by  x  +  (D  —  U)t  =  0.  We  take  this  to  be  the  “zeroth- 
order”  state  upon  which  a  single-mode  transverse  perturbation  is  imposed.  The  position  of 
the  perturbed  shock  can  be  written  as 

xa(t)  =  -(D -U)t  +  Sx(t)eiky ,  (4) 

where  i  =  'J—l,  and  the  real  part  of  the  right  hand  side  of  this  equation  is  implied.  Addition¬ 
ally,  the  wavenumber  k  and  amplitude  5x  satisfy  the  linearity  condition  k8x  <C  1.  Similar 
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Fourier  expansions  hold  for  the  perturbed  hydrodynamic  quantities  behind  the  shock  front. 
Note  that  multi-mode  perturbations  can  be  treated  by  the  inclusion  of  additional  terms  in 
Eq.  (4),  bnt  owing  to  the  linear  independence  of  the  complex  exponential  functions  and 
the  linear  approximation  made  here,  one  can  show  that  doing  so  does  not  affect  our  final 
results.  For  simplicity,  we  therefore  restrict  onr  discussion  to  a  single  wavenumber  only. 
Also  note  that  we  do  not  assume  that  the  amplitude  8x(t)  in  Eq.  (4)  has  the  normal-mode 
form  exp(icat),  with  the  complex  frequency  c o  satisfying  a  dispersion  relation.  While  it  has 
been  shown  that  a  normal-mode  analysis  can  be  used  to  determine  accurately  the  limits  for 
shock  instability  [11,  13],  it  is  also  a  fairly  arduous  and  complicated  approach  to  finding 
a  solution  to  this  problem  [19].  Moreover,  the  normal-mode  formulation  fails  to  capture 
the  correct  temporal  evolution  in  the  stable  regime,  since  the  decay  of  perturbations  on  the 
shock  front  is  not  simply  a  damped  exponential  function  of  time  there.  Our  strategy  for 
analyzing  a  perturbed  isolated  shock  follows  a  different  path.  Namely,  we  seek  a  general 
solution  for  8x(t)  —  or  more  accurately  its  Laplace  transform  —  that  is  valid  under  all 
circumstances.  Such  a  function  contains  complete  information  about  the  evolution  of  the 
rippled  shock  front,  and  once  derived,  can  be  readily  examined  to  determine  the  criteria  for 
unstable  behavior. 

Before  proceeding  further  with  the  summary  of  our  earlier  theory,  it  is  convenient  at  this 
stage  to  introduce  several  dimensionless  expressions.  The  first  is  the  non-dimensional  ripple 
amplitude  g(r),  which  we  write  as 


g(r)  =  k(rj  -  l)Sx/rjr 


(5) 


where  r  =  Dkt/rj  is  a  dimensionless  time.  Next,  we  define  the  additional  quantities 

1  -  Ml 


a2  = 


P  = 

r  = 


Ml  ’ 

1  -h 

2  Mv  ’ 

(1  +  h)rj 

2  Mi  ' 


(6) 

(7) 

(8) 


Note  that  in  terms  of  these  parameters,  the  criteria  h  <  —1  and  h  >  l  +  2Mi  can  be  written 
as  T  <  0  and  (3  <  —  1,  respectively.  Moreover,  the  condition  h  >  hc  is  tantamount  to 
T  >  a2(3.  A  summary  of  the  instability  limits  for  isolated,  planar,  two-dimensional  shock 
waves  as  functions  of  (3,  h,  and  T  is  depicted  graphically  in  Fig.  1.  Throughout  the  remainder 
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of  this  paper,  our  preference  will  be  to  formulate  the  perturbed  shock  problem  in  terms  of 
a,  (3,  and  T,  instead  of  the  more  traditional  set  of  parameters  rj,  Mi,  and  h. 

As  discussed  in  detail  in  Ref.  [27],  one  can  solve  for  the  ripple  amplitude  g{r)  by  adopting 
the  following  procedure.  First,  the  perturbed  hydrodynamic  variables  are  substituted  into 
the  Euler  equations  and  Rankine-Hugoniot  relations  and  linearized.  Then,  using  a  technique 
based  on  Riemann  invariants,  the  resulting  partial-differential  equation  can  be  integrated  in 
time  to  yield  a  Volterra  equation  of  the  second  kind  for  g(r).  The  result  is 


9(t)  =  F(t)+  /  g(0)K(T  -  e)de , 


(9) 


The  kernel  K(q)  appearing  in  the  expression  above  is  given  by 


K(q)  = 


a/3 

1  +  (3 


Ji(aq)  - 


r*aq 


a2/3 


Jo(w)dw 


(10) 


where  Jo  and  J\  are  Bessel  functions  of  order  zero  and  one,  respectively.  For  the  case  of  a 
planar  shock  initially  deformed  into  a  sinusoidal  shape,  the  function  F(t)  in  Eq.  (9)  assumes 
a  particularly  simple  form: 


F(T)  =  9(0)l±^il.  <„) 

Here,  the  symbol  g( 0)  =  F( 0)  denotes  the  initial  (dimensionless)  perturbation  amplitude. 
Taking  the  Laplace  transform  of  Eq.  (9)  leads  to  the  expression 


9l(s)  =  Vs2  +  a2  +  ps  n9^ 

g(0)  s\/s2  +  a2  +  /3s2  +  T 

where  s  is  the  standard  Laplace  variable,  and  the  notation  ^(s)  is  a  convenient  shorthand 
notation  for  denoting  the  transform  of  the  ripple  amplitude.  As  shown  in  Ref.  [27],  this 
function  can  be  manipulated  into  a  sum  of  pairs  of  multiplicative  terms  —  each  of  which  is 
recognizable  as  the  Laplace  transform  of  a  known  function.  LTsing  the  convolution  theorem, 
each  product  then  can  be  inverse  transformed  to  the  time  domain  to  yield  a  final  expression 
for  g(r). 

Following  this  procedure,  it  was  also  demonstrated  previously  that  different  families  of 
stable  solutions  exist,  with  the  envelope  of  oscillations  decaying  asymptotically  in  time  as 
r~3/2.  Although  satisfaction  of  the  inequalities  0  <  T  <  a2 (3  was  tacitly  assumed  in  order 
to  reach  that  conclusion,  it  is  important  to  note  that  such  a  restriction  is  not  a  requirement 
of  our  earlier  work.  Indeed,  the  theory  appearing  in  Ref.  [27]  permits  an  examination  of  the 
nature  of  the  solution  for  any  set  of  parameters  a ,  f3,  and  T,  and  in  this  way,  provides  a 
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means  for  probing  the  stability  of  isolated  planar  shocks.  An  immediately  apparent  approach 
for  conducting  such  an  analysis  is  through  an  inspection  of  the  poles  of  Eq.  (12)  —  a  topic 
that  is  the  subject  of  the  next  section. 


III.  POLES  OF  gL(s)  AND  INSTABILITY  LIMITS 


It  is  well-known  from  control  theory  [28-30]  that  the  poles  of  a  Laplace  transform  de¬ 
termine  the  stability  of  the  time-dependent  system  from  which  it  is  derived.  This  is  a 
remarkable  theorem  of  mathematical  physics,  for  it  implies  that  one  does  not  need  to  com¬ 
pletely  solve  a  dynamical  problem  in  order  to  ascertain  its  late-time  behavior.  Instability  is 
associated  with  poles  that  lie  in  the  right  half  of  the  complex  plane  (i.e.,  have  positive  real 
parts),  while  purely  imaginary  values  imply  stationary  perturbations.  Although  overlooked 
at  the  time  Ref.  [27]  was  published,  the  stability  criteria  appearing  in  Eqs.  (2)  and  (3)  can  be 
readily  inferred  through  the  application  of  this  fundamental  principle  to  Eq.  (12).  Here,  we 
demonstrate  how  this  may  be  accomplished  using  a  simple  method  of  graphical  inspection. 
Additionally,  the  validity  of  the  shock  instability  criteria  is  verified  in  this  section  by  solving 
numerically  the  Volterra  equation  in  Eq.  (9)  for  several  examples  with  different  values  of  a, 
(3,  and  T. 

The  poles  of  Eq.  (12)  occur  where  its  denominator  vanishes.  Thus,  to  investigate  the 
stability  of  g(r),  we  must  solve 

(s  —  ia)1,/2(s  +  ia)1/2  +  /3s  +  T/s  =  0  ,  (13) 


which  has  branch  points  at  s  =  ±ia.  There  are  four  possible  roots  to  this  equation,  which 


In  general,  not  all  of  the  terms  above  will  be  viable  solutions  to  Eq.  (13)  for  a  given  set 
of  shock  parameters  a,  /?,  and  T,  and  one  must  always  confirm  the  validity  of  a  particular 
root  through  explicit  substitution.  For  example,  one  can  show  that  no  solutions  to  Eq.  (13) 
exist  when  the  stability  criteria  0  <  T  <  a2 (3  are  satisfied.  This  is,  of  course,  consistent 


with  the  fact  that  the  damping  of  perturbations  in  the  stable  regime  is  not  an  exponential 
function  of  time  [27].  This  subject  will  be  revisited  in  our  presentation  of  numerical  solutions 
to  the  Volterra  equation  that  follows.  For  the  moment,  however,  we  turn  our  attention  to 
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deriving  analytically  the  conditions  for  which  finite  roots  of  Eq.  (13)  do  exist.  There  are  two 
cases  to  consider,  and  we  discuss  each  separately.  They  are:  (A)  a  root  with  a  positive  real 
part,  which  implies  an  exponentially-growing  mode;  and  (B)  a  root  that  is  purely  imaginary, 
which  corresponds  to  a  mode  that  is  stationary. 

A.  exponentially-growing  modes 

Given  the  particular  form  of  Eq.  (13),  one  can  show  that  a  root  with  a  positive  real  part 
exists  if  and  only  if  the  root  itself  is  purely  real.  We  can  then  look  for  solutions  of  Eq.  (13) 
by  plotting  (s2  +  a2)1/2  and  —(3s  —  T/s  individually  as  functions  of  the  real  variable  s,  and 
determining  under  what  circumstances  these  two  expressions  intersect  for  s  >  0.  Such  a  plot 
is  shown  in  Fig.  2.  The  thick  solid  line  in  this  figure  denotes  the  function  (s2  -t-a:2)1/2,  which 
asymptotes  to  the  linear  function  s  (the  straight  dashed  line  in  Fig.  2)  as  s  — »  oo.  The  thin 
solid  lines  here  represent  the  family  of  curves  given  by  —(3s  —  T/s  for  different  values  of  the 
parameters  (3  and  T.  We  see  that  points  of  intersection  occur  in  two  cases:  (i)  when  T  <  0 
and  (3  >  0,  and  (ii)  when  T  >  0  and  (3  <  —  1.  Both  situations  imply  a  positive  real  root 
of  Eq.  (13),  and  thus  exponentially  growing  (i.e.,  unstable)  solutions.  These  inequalities 
correspond  exactly  to  D’yakov’s  limits  for  the  absolute  instability  of  isolated  planar  shock 
waves  stated  in  Eq.  (2). 

B.  stationary  modes 

Next,  we  wish  to  determine  the  conditions  for  stationary  perturbations  to  exist,  which 
requires  that  we  look  for  purely  imaginary  roots  of  Eq.  (13).  To  do  this,  it  is  useful  to  make 
the  substitution  s  =  ±is,  where  the  variable  s  is  purely  real.  Equation  (13)  then  becomes 

(s2  —  a2)1^  =  —  (3s  +  T/s  .  (14) 

Note  that  in  arriving  at  this  expression,  we  have  used  the  relation  (zLia—ib)1^2  (zLia+ib)1^2  = 
±i(a2  —  62)1/2,  where  a  and  b  are  arbitrary,  but  real  variables.  The  left  and  right  sides  of 
Eq.  (14)  are  plotted  in  Fig.  3.  The  thick  solid  line  denotes  the  function  (s2  —  a2)1/2,  whose 
value  is  not  real  for  values  of  s  less  than  a,  and  asymptotically  approaches  the  linear  function 
s  (the  straight  dashed  line  in  Fig.  3)  as  s  — >  oo.  The  thin  solid  lines  in  this  figure  represent 


the  family  of  curves  given  by  —/3  s  +  T/s  for  different  values  of  the  parameters  /3  and  T.  A 
point  of  intersection  of  these  two  functions  can  occur  only  under  a  particular  circumstance: 
namely,  when  /3  >  —  1  and  —/3a  +  Y/a  >  0,  or  equivalently  Y/a2  >  (3  >  —1.  Consulting 
Fig.  1,  we  see  that  these  inequalities  correspond  exactly  to  the  well-known  criteria  for  the  DK 
instability.  Such  a  solution  implies  the  existence  of  a  purely  imaginary  pole  of  Eq.  (12),  and 
thus  oscillatory,  undamped  solutions.  We  should  remark  that  the  same  criterion  has  been 
derived  previously  by  Wonchuk  and  Cavada  [20]  as  part  of  their  investigation  of  the  unstable 
eigenmode  spectrum  of  shocks  launched  from  rippled  pistons.  In  that  study,  the  authors 
employed  a  graphical  technique  similar  to  the  one  adopted  here  to  analyze  the  stability 
properties  of  Zaidel’s  Bessel-series  expansion  for  the  downstream  pressure  held  [31]. 

C.  numerical  solutions  of  the  Volterra  equation 

Figure  4  shows  five  examples  of  the  evolution  of  the  amplitude  of  a  linear  sinusoidal 
perturbation  on  a  shock  front  obtained  by  numerically  solving  the  Volterra  equation  in 
Eq.  (9)  for  different  values  of  a2,  /3,  and  Y  [32],  In  performing  these  calculations,  one 
must  ensure  that  the  algorithm  employed  to  compute  the  values  of  Bessel  functions  of  large 
argument  does  so  with  sufficient  accuracy.  The  curves  in  Fig.  4  are  labeled  “a”  through  “e,” 
and  the  parameters  for  each  case  are  listed  in  Table  I,  along  with  the  corresponding  values 
of  the  poles  of  Eq.  (12).  Note  that  the  curves  labeled  “a”  and  “b”  satisfy  the  inequalities 
/3  <  —  1  and  Y  <  0,  respectively,  and  are  clearly  unstable  in  an  absolute  sense.  Accordingly, 
both  have  a  pole  that  is  real  and  positive.  Each  of  the  curves  labeled  “c”  and  “d”  satisfies 
Y/a 2  >  /3  >  —1,  and  each  exhibits  a  stationary  instability  with  perturbations  that  persist 
indefinitely;  also  note  that  both  of  these  curves  have  purely  imaginary  poles,  as  well  as  a 
negative  real  one,  corresponding  to  an  exponentially-damped  transient  solution.  Curve  “e” 
is  an  example  of  a  stable  case;  it  has  no  poles  and  a  perturbation  amplitude  that  damps 
asymptotically  as  r~ 3/2. 

IV.  CONCLUSIONS 

The  present  work  is  an  outgrowth  of  a  previous  theoretical  investigation  [27]  on  the 
effect  of  small  disturbances  to  an  isolated,  planar,  two-dimensional  shock  wave  propagating 
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through  a  fluid  with  an  arbitrary  EOS.  In  that  study,  an  explicit  expression  for  the  time- 
dependent  Fourier  coefficient  associated  with  a  linear  single-mode  perturbation  on  the  front 
was  obtained  by  deriving  an  equation  for  the  Laplace  transform  of  the  ripple  amplitude. 
Although  overlooked  at  the  time  Ref.  [27]  was  published,  the  theory  presented  there  is 
sufficiently  general  to  permit  an  examination  of  the  nature  of  the  solution  for  any  set  of 
parameters  a,  /3,  and  T,  and  thus  provides  a  useful  method  for  analyzing  the  stability  limits 
of  isolated,  planar,  two-dimensional  shock  waves.  As  we  have  seen,  a  particularly  elegant 
approach  for  conducting  such  an  investigation  is  through  an  inspection  of  the  poles  of  the 
Laplace-transformed  ripple  amplitude  —  an  expression  for  which  appears  in  Eq.  (12).  By 
analyzing  the  poles  of  this  Laplace  transform,  we  have  demonstrated  that  two  instability 
classes  exist  for  isolated  planar  shocks:  one  in  which  perturbations  grow  exponentially  in 
time,  and  the  other  in  which  disturbances  are  stationary.  The  necessary  criteria  for  the  two 
classes  of  instabilities  to  occur  can  be  stated  as  T  <  0  or  f3  <  —1,  and  Y/a2  >  (3  >  — 1, 
respectively.  These  results  agree  with  those  derived  by  D’yakov  [11]  and  Kontorovich  [16] 
by  more  laborious  means,  and  serve  as  an  important  addendum  to  our  earlier  work.  In 
addition,  we  have  shown  in  this  paper  that  such  limits  are  consistent  with  several  examples 
obtained  by  numerically  solving  the  Volterra  equation  in  Eq.  (9)  for  different  combinations 
of  the  constants  cc,  /3,  and  T. 

We  wish  to  emphasize  that  the  instability  limits  derived  in  this  study  apply  to  a  particular 
category  of  shock  waves.  This  is  a  consequence  of  the  fact  that  in  arriving  at  those  criteria, 
several  approximations  were  made,  two  of  which  warrant  additional  comment  here.  The 
first  is  the  assumption  that  the  linear  sinusoidal  deformation  to  the  otherwise-planar  shock 
front  occurs  suddenly  at  r  =  0,  and  we  make  no  inquiry  as  to  the  actual  physical  mechanism 
( e.g .,  a  constricted  boundary  channel  [18])  by  which  such  a  perturbation  might  arise.  As  a 
result,  the  evolution  of  the  rippled  shock  can  be  treated  in  the  context  of  an  initial- value 
problem  with  periodic  transverse  boundary  conditions.  A  second  simplifying  assumption  is 
that  the  shock  is  “isolated,”  i.e.,  it  is  far  enough  away  from  the  “piston”  driving  it  that  the 
transit  time  of  a  sound  wave  between  these  two  surfaces  is  much  longer  than  the  time  of 
interest  for  this  problem.  Such  an  arrangement  is,  of  course,  an  idealization  and  one  may 
ask  how  the  instability  limits  would  change  for  a  more  realistic  scenario  in  which  acoustic 
waves  impinge  upon  —  and  reflect  off  —  a  shock  front  from  the  compressed  fluid  region.  To 
a  large  extent,  this  remains  an  open  question.  We  should  remark,  however,  that  in  the  case 
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of  the  DK  instability,  Wouchuk  and  Cavada  [20]  have  shown  that  additional  eigenmodes 
of  oscillation  are  excited  on  the  front  as  the  result  of  a  reflecting  surface  behind  it.  There 
is  also  evidence  to  suggest  that  the  absolute  instability  limits  appearing  in  Eq.  (2)  require 
modification  in  the  presence  of  such  a  surface  [12,  33]. 
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Case 

P 

r 

poles 

comments 

a 

-1.100 

+2.553 

+6.216,  -0.896 

positive  real  pole;  exponentially  growing  solution  (/?  <  —1) 

b 

+1.064 

-0.100 

+0.059,  -4.669 

positive  real  pole;  exponentially  growing  solution  (T  <  0) 

c 

-0.500 

+2.553 

-1.040,  ±2.836 i 

purely  imaginary  poles;  stationary  perturbation  (T/a2  >  /3  >  —1) 

d 

+0.500 

+2.553 

-1.695,  ±1.739 i 

purely  imaginary  poles;  stationary  perturbation  (T/a2  >  (3  >  —1) 

e 

+1.064 

+2.553 

— 

no  poles;  stable  solution  (perturbation  damps  as  r~3/2  as  r  — ■>  oo) 

TABLE  I:  Parameters  for  the  five  curves  shown  in  Fig.  4.  The  value  of  a2  in  all  cases  is  2.667. 
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LIST  OF  FIGURE  CAPTIONS 


FIG.  1:  Unstable  regions  of  isolated,  planar,  two-dimensional  shock  waves  in  terms 
of  various  dimensionless  parameters.  Definitions  are  given  in  the  text.  Note  that  the 
parameters  (3  and  Y  are  not  independent  of  each  other:  if  Y  <0,  then  (3  >  0,  and  if 
(3  <  —1,  then  Y  >  0.  For  Y  >  0,  however,  the  parameter  (3  may  be  either  positive  or  negative. 

FIG.  2:  Plots  of  (s2  +  cc2)1/2  (thick  line)  and  the  family  of  curves  given  by  —/3s  —  Y/s  (thin 
solid  lines  with  different  values  of  [3  and  T)  as  a  function  of  the  real  variable  s.  Points  of 
intersection  on  this  graph  denote  the  existence  of  positive  real  poles  of  Eq.  (12),  and  thus 
exponential  growing  (i.e.,  unstable)  solutions.  This  can  occur  in  two  cases:  (i)  when  Y  <  0, 
and  concomitantly  (3  >  0,  and  (ii)  when  (3  <  —1,  and  concomitantly  Y  >  0. 

FIG.  3:  Plots  of  (s 2  —  a2)1/2  (thick  line)  and  the  family  of  curves  given  by  —/3s  +  T/s  (thin 
solid  lines  with  different  values  of  (3  and  T)  as  a  function  of  the  real  variable  s.  Points  of 
intersection  on  this  graph  denote  the  existence  of  purely  imaginary  poles  of  Eq.  (12),  and 
thus  oscillatory,  undamped  solutions.  This  can  occur  when  (3  >  —  1  and  the  inequality 
—(3a  +  Y /a  >  0  is  satisfied.  Note  that  the  latter  condition  can  be  rewritten  as  Y /a2  >  (3. 

FIG.  4:  Five  examples  showing  the  evolution  of  the  amplitude  of  a  linear  sinusoidal 
perturbation  on  a  shock  front  obtained  by  numerically  solving  the  Volterra  equation  in 
Eq.  (9).  Each  case  corresponds  to  a  different  set  of  dimensionless  parameters  a2,  (3,  and  T; 
see  Table  I. 
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